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Abstract 

This technical note describes the application of saddle-point integra¬ 
tion to the asymptotic Fourier analysis of the well-known Coo “bump” func¬ 
tion e~^~ x ■* l(_i,i), deriving both the asymptotic decay rate k~ 3 C e ~Ck 

of the Fourier transform F(k ) and the exact coefficient. The result is 
checked against brute-force numerical integration and is extended to gen¬ 
eralizations of this bump function. 


1 Background 


In our work involving optimization of slow “taper” transitions for waveguide 
couplers and also for our work on exponential absorbing boundary layers in 
wave equations 00 - we encountered a problem that could be expressed in 
terms of the Fourier transform of smooth functions with compact support. In 
particular, we needed to analyze the asymptotic rate of decay of the Fourier 
coefficients. This was straightforward for functions with simple discontinuities 
(via integration by parts HHU) or with sim P le P° les ( via colour integration (5j 
©)■ However, it turned out that some of the most interesting functions were Coo 
(infinitely differentiable) functions with compact support, and these functions 
have essential singularities that resist those methods. In this note, we describe 
how the asymptotic Fourier transforms of such functions can be analyzed with 
the help of saddle-point integration 17 Q A similar approach has been applied to 
various functions in the related context of Chebyshev polynomial series (8 10. 

In particular, we will look at Coo functions on R. with compact support 
[—1,1]. The canonical example of such a function is the symmetric “bump:” 


/(d = 


a; €(-1,1) 
otherwise 


The actual functions we are interested may be more complicated than this [e.g., 
f(x) multiplied by some analytic function], but their analysis is similar to that 

1 We are indebted to our colleagues, Prof. Martin Bazant and especially Prof. Hung Cheng 
at MIT, for their helpful suggestions in this matter. 
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of f(x): the key point is that we have essential singularities at x = — 1 and 
x = +1, and these singularities determine the asymptotic behavior of the Fourier 
transform and similar integrals. 

We should also note that the space of functions with compact support 
plays an important role in the theory of generalized functions (distributions), 
where they are the “test functions” that are the domain of the distributions. 
In this context, it has been proven that the Fourier transform of any such 
test function is an entire function (analytic over the whole complex plane) and 
diverges at most exponentially fast off the real axis [II]. From the fact that 
the functions are infinitely differentiable, it also immediately follows that their 
Fourier transforms go to zero along the real axis faster than the inverse of any 
polynomial [3|4|. Here, however, we want to know precisely how fast the Fourier 
transform decays, and with what coefficient. 


2 Saddle-point Fourier integration of fix) 


We wish to compute the asymptotic behavior of the Fourier transform 


F(k) 


f{x)e ikx dx = 2 Re 


ikx— 

e 


1 



Jo 


for Re k 1. (Without loss of generality, we can restrict ourselves to real 
k > 0.) To do this, we will employ a saddle-point integration: we will look for the 
x at which the exponent is stationary, and expand the exponent approximately 
around this point. For large k, this saddle point of the stationary exponent 
will dominate the integral, and this will give us the asymptotic behavior. It 
will turn out that the saddle point occurs for a complex x, however, so we will 
need to deform the integration contour within the complex plane to exploit this 
approach. 

It is clear (and this will also be justified a posteriori ) that the biggest con¬ 
tributions must come near the singular point x = 1. Exactly at the singular 
point, however, the integrand is zero, so (perhaps) contrary to our intuition the 
endpoints per se are not important. Because we are worrying about points near 
the endpoint, however, it is convenient to change variables t = 1 — x and write 

F(k)= 2Re / e 9W dt 
Jo 

with 

q(t) = ik — ikt — ~ ik — ikt - V O(f), 

yw (2-t)t a 4 w 

where the last approximation is for |t| <C 1, which is valid (it turns out) at our 
saddle point. The saddle point is the t = t 0 where g'(t) = 0 = —ik + l/2t 2 , 
which by inspection i^] to = \J\j2%k. Note that for large i; > 1 we obtain 

1 g' also vanishes at t = — y l /2ik. but we cannot deform our integration contour to a point 
with negative Ret: for Ret < 0, our integrand ~ e —l/ 2i blows up. 
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|to| < 1, justifying our approimation above. Now, we write 


g(t) ~ g{t o) + a 2 — ^o ) 2 



To actually do this integral, we need to deform our integration contour in the 
complex plane to lie along a line t = u/y/i for real u near u = 0, so that the 
saddle point t = to lies along our integration path (at u = ^/l/2fc)j^] (This 
deformation is not a problem since we don’t have any singularities except at 
t = 0 and t = 2, and the integrand vanishes as t —>■ 0 for Ret > 0.) In this case, 
our integral becomes (approximately) a Gaussian integral, since: 


g{u/\fi) =ik- - 


y/2 ik — \[2ik 3 / 2 




Recall that the integral of J^° e~ au du = \Jn/a as long as Re a > 0, which is 
true here. Note also that, in the limit as k becomes large, the integrand becomes 
zero except close to u = \Jl/2k , so we can neglect the rest of the contour and 
treat the integral over u as going from — oo to oo. (Thankfully, the width of 
the Gaussian Art ~ k ~ 3 / 4 goes to zero faster than the location of the maximum 
uo ~ A; -1 / 2 , so we don’t have to worry about the u = 0 origin.) Also note that 
the change of variables from t to u gives us a \J—i Jacobian factor. Thus, when 
all is said and done, we obtain the exact asymptotic form of the Fourier integral 
for k 1: 


m 


2 Re 


' ~ l7T c ik-^-V2ik 

V2ik 3 / 2 


Since y/2i = 1 + i, this means that the Fourier coefficients decay proportional to 
fc _3 / 4 e _v/ ^, which is consistent with the expected faster-than-polynomial decay. 

To check this result, we can compare the above formula with an exact nu¬ 
merical evaluation of the Fourier integral. Numerical integration was performed 
using adaptive Clenshaw-Curtis quadrature, specialized for a cos (kx) oscillatory 
weight function, from the GNU Scientific Library (adapted from QUADPACK). 
The resulting |i r (A;)| for 0 < k < 150 is plotted in figure [l] As can be seen from 
the figure, the asymptotic approximation is an excellent match for the exact 
result, with errors under 10% for frequencies k as small as 4. 

3 Although t = u/y/i is not the path of steepest descent around to, as would be prescribed 
by the usual saddle-point method (a.k.a. the method of “steepest descent”), it is at least a 
path of descent (the integrand is decaying along that path) [7j. This is sufficient for us to 

apply our Gaussian integral approximation. 
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Figure 1: “Exact” numerical quadrature of |F'(A:)| (red line) and asymptotic 
saddle-point integration (blue line) for the Fourier transform of the bump 
function f(x). (They match so well that the curves are difficult to distinguish.) 
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3 Generalized bumps 

As a further warm-up for the more complicated problems we may want to solve 
later, let’s look at the asymptotic Fourier transform of a generalization of our 
bump function f(x): 


fa, id(x) 


e v-J)*- 1 x €(- 1 , 1 ) 
0 otherwise 


so that f{x) = f' 2 ;i ■ Again, we look at the Fourier transform with changed 
variables x = 1 — t, i.e. f* / a>/ 3(1 — t)e lk ~ lkt dt, and expand the exponent g(t) 
around t = 0 : 


g{t) = ik — ikt 


_ 

(2 - f ) 0 " 1 


ta - 1 


ik — ikt 


0 

(2t) a ~ 1 


+ 0{t 2 ~ a ). 


Our derivative is then g'(t) = 0 « —ik + f3{a — 1)/2 Q 1 t a at 

t 0 = [2 1 _ “/3(a - 1 )/ik] 1,a . 

That is, to ~ k ~ x ! a , and therefore g(to) ~ asymptotically. 

The second derivative is g" (to) = — 2 1 -Q /3a(a— l)/fQ +1 = — i( a + 1 )/ a 2Ak^ a+1 ^ a 
where 

A = a[2(3{a-l)}~ 1/o: . 

Again, we’ll choose a contour t = u/i 1 /“, in which case we find: 


g(t) « <?(<o) - * (a " 1 )/Q AA:(“ +1 )/“(u - u 0 ) 2 , 


which is a path of descent so we can perform a Gaussian integral. The final 
answer for the integral, including the Jacobian factor for dt = du/i 1 ^ 01 , is then 


F{k) « 2 Re 


ik—ikto — 


(ifc)(° +1 )/“A 


[(2 —t 0 )to] Q 


with to and A given above. Notice that for our exponent g(to) we use the exact 
form of g(t) and not its approximation for small t. The reason is that, for 
a > 2 , terms of 0(t 2 ~ a ) do not go to zero, and therefore make a non-negligible 
multiplicative contribution to the amplitude F{k) even though they do not affect 
the saddle-point integration. 

Numerical tests seem to confirm the accuracy of this F(k) formula, although 
for a > 3 we start to have more difficulty with the numerical quadrature. As one 
might have expected, increasing either a or f3 makes F(k) decay more rapidly. 
To perform such comparisons more carefully, we should typically normalize by 
Jo fa, ft(x)dx or similar, to ignore effects due simply to the fact that the integrand 
is getting smaller overall. 
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